Experimental observation of structural crossover in binary mixtures of colloidal hard 

spheres 
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Using confocal-microscopy we investigate the structure of binary mixtures of colloidal hard spheres 
with size ratio q = 0.61. As a function of the packing fraction of the two particle species, we observe 
a marked change of the dominant wavelength in the pair correlation function. This behavior is in 
excellent agreement with a recently predicted structural crossover in such mixtures. In addition, the 
repercussions of structural crossover on the real-space structure of a binary fluid are analyzed. We 
suggest a relation between crossover and the lateral extension of networks containing only equally 
sized particles that are connected by nearest neighbor bonds. This is supported by Monte-Carlo 
simulations which are performed at different packing fractions and size ratios. 

PACS numbers: 82.70.Dd, 61.20.-p 



Most systems in nature and technology are mixtures 
of differently sized particles. Each distinct particle size 
introduces another length scale and its competition gives 
rise to an exceedingly rich phenomenology in compari- 
son with single-component systems. Already the simplest 
conceivable multi-component system, i.e. a binary mix- 
ture of hard spheres, exhibits interesting and complex 
behavior. Just a few examples include entropy driven 
formation of binary crystals [H HI [3], frustrated crys- 
tal growth [4], the Brazil nut effect [5], glass- formation 
^, 7 and entropic selectivity in external fields [8 . Al- 
though interaction potentials in atomic systems are more 
complex than those of hard spheres, the principle of vol- 
ume exclusion is ubiquitous and thus always dominates 
the short-range order in liquids [9]. Accordingly, hard 
spheres form one of the most important and successful 
model systems in describing fundamental properties of 
fluids and solids. It has been demonstrated that many 
of their features can be directly transferred to atomic 
systems where fundamental mechanisms are often ob- 
structed by additional material specific effects [10]. Bi- 
nary hard sphere systems are fully characterized by their 
size ratio q = (Ts/o'b with the diameters of the small 
(S) and big (B) spheres and the small and big sphere 
packing fractions 7^5, r]B^ respectively. 

The pair-correlation functions, gij{r), are the central 
measure of structure in fluids; they describe the probabil- 
ity of finding a particle of size i at distance r from another 
particle of size j. It is well known that all pair-correlation 
functions in any fluid mixture with short-ranged inter- 
actions (not just hard spheres) exhibit the same type 
of asymptotic decay, which can be either purely (mono- 
tonic) exponential or exponentially damped oscillatory 
(pT| and references therein). This prediction, which is 
valid in all dimensions, suggests that all pair-correlation 
functions decay with a common wavelength and decay 
length in the asymptotic limit. For binary hard- sphere 
mixtures where r]B ^ Vs or r]s ^ Vb^ this is rather obvi- 
ous since the system is dominated by either big or small 



particles. The pair-correlation functions will asymptot- 
ically oscillate with a wavelength determined either by 
ctb {vb > Vs) or as {r]s > Vb)- Rather surprising is 
that the above statement is also valid for all other rela- 
tive packing fractions where the system is not dominated 
by particles of a single size ([UlIIl])- Accordingly, in the 
asymptotic limit the {tjs^Vb) phase diagram is divided 
by a sharp crossover line where the decay lengths of the 
contributions to gij (r) with the two wavelengths become 
identical. Below and above this line, however, the pair- 
correlation function is either determined by the diameter 
of the small spheres or that of the big spheres [13 . 

Despite the generic character of structural crossover 
and the close relationship between structural and me- 
chanical properties, this effect has not been observed in 
experiments as the asymptotic limit is difficult to reach 
in scattering experiments on atomic and molecular liq- 
uids. However, recent calculations suggest that struc- 
tural crossover is already detectable at relatively small 
distances [12^. Because colloidal particles are directly 
accessible in real space, such systems provide an oppor- 
tunity to explore the structure of binary fluids and to 
investigate structural crossover experimentally. 

As colloidal suspension we used an aqueous binary 
mixture of small melamin particles {as = 2.9/im) and 
big polystyrene spheres {as = 4.8/im). Addition of 
salt screens residual electrostatic interactions thus lead- 
ing to an effective hard sphere system. Since melamin 
has a higher density {pm = 1.51^/cm^) than polystyrene 
{Pp = l.Obg/cm^) the sedimentation velocities are sim- 
ilar and, hence, we obtain a homogeneous system after 
mixing. The suspension was contained in a cylindrical 
sample cell with a silica bottom plate to allow optical 
imaging with an inverted confocal microscope in reflec- 
tion mode (Leica TCS SP2). From the images, particle 
positions were obtained with digital video microscopy 
p^. Strong layering at the bottom wall allowed us to 
image only the first two-dimensional bottom layer of the 
three-dimensional system. We define the packing fraction 
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Figure 1: Different paths with constant total packing fraction 
Tf — rfs -\- TfB in the {tis^tib)- plane. Experimental data (open 
symbols: 77 = 0.72, q = 0.61) are sorted into ten bins. The 
bin size is indicated by the 'error bars'. Closed symbols cor- 
respond to the MC-simulations (A: rf — 0.62, q — 0.4) and 
(•: Tj = 0.57, q = 0.5). For convenience all samples are la- 
beled with numbers increasing in the direction indicated by 
the arrows. 



as Tji = 7r(j|/4, with pi the number density of component 
i. Variation of the relative packing fractions of the par- 
ticles was achieved by addition of small particles to a 
suspension of big spheres (Figjl]). Thus, the total pack- 
ing fraction in the two-dimensional bottom layer remains 
constant for all samples: 77 = 0.72. In the following we 
will refer to the different samples by the sample numbers 
(No.) as given in Figjl] 

Typical snapshots of the system for different packing 
fractions of big and small particles are shown in Figs|2}^- 
C. The images demonstrate how the structure of the bot- 
tom layer changes from being rich in small particles (No. 
1, Figj2]A.) to being rich in big particles (No. 10, Figj2]C). 
Figj2^ (No. 5) corresponds to about the same number 
density of small and big spheres. In order to analyze the 
samples for a possible structural crossover, we calculated 
the pair correlation function from the determined particle 
positions. To minimize statistical noise we did not distin- 
guish between big and small spheres. This is justified be- 
cause the crossover has been predicted to be visible in all 
pair-correlation functions and thus also in any linear com- 
bination [11| |12^ . The dominating wavelength in the os- 
cillations is identified by computing the total correlation 
function htot{r) = Y^^j XiXjhij{r) = Y,.. XiXj[gij{r) -I], 
with the mole fraction Xi = pi / pi of component i p!2l . 
Figj2|l) exemplarily shows In \ htot{'f^) \ for samples No. 1,5, 
and 9. Note that in this representation the oscillation 
wavelength is halved. The correlation functions of sam- 
ples No.l and 9 clearly oscillate with a single wavelength, 
respectively, given by ^ (Jb/^^ and ^ (Js/2. In contrast, 
sample 5 does not show a dominating wavelength but 
an interference of different length scales which is typical 
near the structural crossover. It is important to mention. 
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Figure 2: A-C) Typical snapshots of the bottom layer of a 
binary mixture observed with a confocal microscope used in 
reflection mode. The mixtures correspond to sample 10 (A), 5 
(B) and 1 (C). The field of view is 40x40/im^. D) Logarithmic 
plot of the total correlation functions htot(r) for the experi- 
mental binary mixtures with 77 = 0.72 zb 0.04. Correlation 
functions are plotted for sample numbers 1,5 and 9 (compare 
Figjl]) and are shifted in vertical direction for clarity. The hor- 
izontal bars correspond to (7^/2 and crs/2, respectively. E) 
Fourier-transforms of htot{r) for the experimental data points 
(compare Fig. 1). Vertical lines indicate the wave vectors k 
corresponding to the diameters of the small (S) and big par- 
ticles (B), respectively, (color online). 



that this intermediate behavior is only observed for sam- 
ples No. 5 and 6, i.e. only for about 10% of the entire 
range over which tjb and rjs was varied. The experimen- 
tally identified crossover-region is in excellent agreement 
with the theoretically calculated value of 775 « 0.3 at 
those size ratios, which were determined from the decay 
of the pair correlation functions calculated within density 
functional theory in the test particle limit [T5^. Fig.2E 



Figure 3: Visualization of the different bond- types as determined by a Delaunay triangulation: big-big (black), big-small 
(yellow) and small-small (red). Different plots correspond to the sample numbers as indicated in Figjl] The field of view is 
180 X 180/im^ 



shows the Fourier transforms of htot{^) for all samples 
where the rather sudden change of the dominating wave- 
length is seen more clearly [16]. At small and high pack- 
ing fractions, the correlations are clearly dominated by 
frequencies corresponding to either small or large parti- 
cles (vertical lines) while around sample No. 5 hardly any 
dominating frequency is observed. This experimentally 
confirms structural crossover as well as its occurrence at 
finite particle distances. 

So far, structural crossover has been discussed in terms 
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Figure 4: Averaged radii of gyration {RD (normalized to L/2 
with the size of the field of view) of networks formed by 
large (solid symbols) and small particles (open symbols) as a 
function of the sample number for A) the experimental data, 
B) the MC-simulations at 77 = 0.57 and q — 0.5 and, C) the 
MC-simulations at 77 = 0.62 and q — 0.4. The correspond- 
ing packing fraction of small particles rjs is indicated as well. 
The grey area and the dashed line respectively indicate the 
crossover as inferred from the correlation functions and from 
density functional theory, (color online). 



of pair correlation functions, i.e. spatially averaged quan- 
tities. Since our experiments naturally provide detailed 
structural information, we investigate what the reper- 
cussions are of the structural crossover on the real-space 
structure. We first subjected a Delaunay triangulation to 
the set of particle centers and identified nearest-neighbor 
bonds between big-big (black), big-small (yellow), and 
small-small (red) particles, respectively (see Figlsl. As 
observed in Fig|3j sample 1 predominantly consists of big- 
big bonds which form a large network spreading across 
the entire field of view. With increasing sample No., i.e. 
increasing 775, the number of small-small bonds increases, 
which leads to fragmentation of the big-big network into 
smaller, randomly distributed patches. At large sample 
numbers, the role of big and small particles is inverted 
and small-small bonds form a network spanning the en- 
tire area (No. 10). Having distinguished between differ- 
ent bond- types, a natural and well-known measure of the 
spatial extend of a network formed by particles of size 
i at positions {k = 1 . . .rf) \s given by the radius of 

gyration R'g = with R'q the cen- 

troid position of the network. Computing this quantity 
for all, say N^^ networks formed by connected particles 
of size finally yields a weighted averaged radius of gyra- 
tion {Rg) = ^m=i ni{'m)Rg{m) where N'^ denote the 
total number of particles i. We calculated (Rg) for net- 
works consisting of connected big or small particles and 
plotted these values for our experimental data in Figj4j\ 
as a function of the sample number. At small and high 
sample numbers the quantities saturate while a relatively 
sharp transition with an intersection point occurs around 
sample 6. This location is indeed in very good agreement 
with the crossover transition as determined from the cor- 
relation functions in Fig j2] and density functional theory 
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(also indicated in Fig. |4j^). This suggests that the struc- 
tural crossover corresponds to a competition between the 
sizes of networks consisting of connected big or small par- 
ticles, respectively. 

As structural crossover is also predicted for other size 
ratios and packing fractions, we use Monte-Carlo (MC) 
simulations to test our findings for more dilute systems 
with size ratios q = 0.5 and q = 0.4. The corresponding 
paths through the phase diagram (see closed symbols in 
Figjl]) were obtained from 2-dimensional simulations with 
a fixed number of particles of about < N < 3000 for 
both species and box areas of about 1500cr^ employing 
periodic boundary conditions. From the configurational 
snapshots we first determined the region of crossover by 
analyzing htoti'^) (the correlation functions are sampled 
using 10^ MC cycles per particle). Then, we performed 
the above described Delaunay triangulation to calculate 
(Rg) for networks of connected big or small particles, re- 
spectively. The corresponding radii of gyration are plot- 
ted in Figj4^ and C and show a similar behavior as in the 
experiment. Again, the intersection points are consistent 
with the crossover region as inferred from the correlation 
functions and DFT calculations. Note that the crossover 
region sensitively depends on the size ratio and packing 
fractions. Both the experiment and Monte- Carlo sim- 
ulations show that structural crossover is accompanied 
by a pronounced change in the typical size of networks 
consisting of connected big and small particles. By in- 
troducing small particles into a system of big spheres. 



connections between big particles are broken and, at the 
same time, connections between small particles are made. 
This sensitively affects the typical size of networks con- 
taining connected, equally-sized particles and thereby the 
chance of finding another particle with the same size at a 
relatively large distance. Consequently, the change from 
(Rg) > {R^) to {R^) > {Rg) (and vice versa) provides 
a simple real-space argument why the oscillation wave- 
length of the gij(r) in the asymptotic limit is either set 
by gb or gs- 

We have experimentally demonstrated the structural 
crossover in a binary colloidal hard sphere system. Fur- 
thermore, we show that structural crossover is strongly 
coupled to the size of networks containing connected 
equally-sized particles only. Going across the structural 
crossover, the size ratio of such networks comprised by 
either connected big or small particles is reversed. We be- 
lieve this real-space configurational picture of structural 
crossover is not just applicable to binary hard spheres, 
as structural crossover is a generic feature of mixtures 
with competing length scales. Moreover, it shows inter- 
esting similarities with force chains in granular matter 
[17] and glassy systems [H [7[ |18] of dissimilar sized parti- 
cles. Therefore, our finding may help to gain more insight 
into structure-related properties in binary systems at an 
universal level. 
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